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PREFACE 


-  This  r«port  describes  the  work  completed  during  the  first  year  of  a 
three-year  investigation  into  the  feasibility  of  using  in-situ 
.  observations  of  the  ionosphere  from  the  DMSP  SSIES  sensors  to 
calculate  parameters  which  characterize  ionospheric  scintillation 
effects.  This  work  is  part  of  a  larger  effort  with  an  overall 
objective  of  providing  the  USAF  Air  Weather  Service  with  the 
capability  of  observing  ionospheric  scintillations,  and  the  plasma 
density  irregularities  which  cause  the  scintillations,  in  near  real¬ 
time  and  updating  models  of  ionospheric  scintillation  with  these 
observations. 


1 .  INTRODUCTION 


Many  modern  military  systems  used  for  communications,  command  and 
control,  navigation,  and  surveillance  depend  on  reliable  and 
relatively  noise-free  transmission  of  radiowave  signals  through  the 
earth's  ionosphere.  Small-scale  irregularities  in  the  ionospheric 
density  can  cause  severe  distortion,  known  as  radiowave  scintillation, 
of  both  the  amplitude  and  phase  of  these  signals,  A  basic  tool  used 
in  estimating  these  effects  on  systems  is  a  computer  program,  WBMOD, 
based  on  a  single-scatter  phase-screen  propagation  model  and  a  number 
of  empirical  models  of  the  global  morphology  of  ionospheric  density 
irregularities.  An  inherent  weakness  of  WBMOD  is  that  the 
Irregularity  models  provide  median  estimates  for  parameters  with  large 
dynamic  ranges,  which  can  lead  to  large  under-  and  over-estimation  of 
the  effects  of  the  ionospheric  irregularities  on  a  system. 

One  solution  to  this  problem,  at  least  for  near  real-time 
estimates,  is  to  update  the  WBMOD  irregularity  models  with 
observations  of  the  various  parameters  modeled.  One  proposed  source 
for  these  observations  is  from  the  in  situ  plasma  density  monitor  to 
be  flown  on  the  Defense  Meteorology  Satellite  Program  (DMSP) 
satellites.  This  study  is  designed  to  assess  the  applicability  of 
this  data  set  to  real-time  updates  of  the  WBMOD  models.  There  are  two 
primary  objectives: 

(1)  Develop  and  refine  techniques  for  generating  estimates  of 
parameters  which  characterize  ionospheric  scintillation  from  in  situ 
observations  of  the  ionospheric  plasma  from  the  DMSP  SSIES  sensors. 

(2)  Determine  if  the  parameters  calculated  from  the  SSIES  data 
can  be  used  to  determine  the  scintillation  effects  on  a 
transionospherlc  radiowave  signal. 

This  report  describes  the  results  obtained  during  the  first,  year 
of  the  study.  The  focus  during  this  year  was  on  developing  the 
techniques  to  be  used  for  calculating  scintillation  parameters  from 
the  SSIES  data  and  making  determinations  of  the  uncertainties  involved 
in  these  calculations. 


2 .  BACKGROUND 


The  propagation  model  used  in  the  WBMOD  program  (based  on  weak- 
scatter  phase-screen  theory!1! )  characterizes  the  ionospheric  electron 
density  irregularities  which  cause  scintillation  via  eight  independent 
parameters!2!: 

(1)  The  irregularity  axial  ratio  along  the  direction  of  the 
ambient  geomagnetic  field,  a. 

(2)  The  irregularity  axial  ratio  perpendicular  to  the  direction 
of  the  ambient  geomagnetic  field,  b. 

(3)  The  angle  between  sheet-like  irregularity  structures  and 
geomagnetic  L  shells,  f>. 

(4)  The  height  of  the  equivalent  phase  screen  above  the  earth's 
surface,  hp. 

(5)  The  in  situ  irreqularity  drift  velocity,  v^. 

(6)  The  outer  scale  of  the  irregularity  spectrum,  iv0. 

(7)  The  slope  of  a  power-law  distribution  which  describes  the 
one-dimensional  power  density  spectrum  ( PDS )  of  the  irregularities,  q. 

(8)  The  height-integrated  strength  parameter,  C^L. 

The  first  three  parameters  (a,  b,  and  6)  and  the  direction  of  the 
ambient  geomagnetic  field  specify  the  propagation  geometry,  while  the 
last  three  (iv0,  q,  and  C^L)  specify  the  spectral  characteristics  of 
the  Irregularities. 

It  may  be  possible  to  obtain  astimates  for  the  values  of  three  of 
these  parameters  from  the  DMSP  SSIES  sensors:  (from  the  SSIES  ion 

Drift  Meter  (DM)),  and  q  and  C^L  (from  the  SSIES  ion  Scintillation 
Meter  (5M)).  In  this  study,  we  will  focus  on  the  estimation  of  C^L 
from  this  data  set  and  consider  q  and  y^  only  in  terms  of  the  effects 
of  uncertainties  in  these  parameters  on  the  estimates  of  C^L.  Of  the 
eight  parameters,  C^L  varies  the  most  as  a  function  of  location  and 
time,  and  has  the  most  profound  effect  of  the  accuracy  of  estimates  ot 
scintillation  levels  made  by  the  WBMOD  model. 

In  the  phase-screen  propagation  theory  used  in  WBMOD, Is!  the  C^I. 
parameter  is  actually  the  product  of  two  parameters:  the  three 


dimensional  spectral  "strength"  of  the  electron  density  irregularities 
at  a  scale  size  of  1000km  (related  to  the  structure  constant  used  in 
classical  turbulence  theory);  and  L,  the  thickness  of  the  irregularity 
layer'*  The  models  in  WBMOD  ware  obtained  from  analysis  of  phase 
scintillation  data  from  the  WIDEBAND  and  HiLat  satellites,  which  will 
provide  estimates  of  the  height-integrated  value  of  C^L  rather  than 
independent  measures  of  and  L.  Because  of  this,  the  model  was 
developed  for  C^L  rather  than  for  and  L  separately. 

The  calculation  of  an  estimate  of  the  C^L  parameter  from  topside 
in  situ  ion  density  observations  requires  two  operations.  First,  an 
estimate  of  at  the  satellite  altitude  is  made  from  a  finite-length 
time  series  of  density  measurements.  Second,  the  estimate  of  is 
converted  to  an  estimate  of  C^L  in  some  fashion  which  will  account  for 
both  the  thickness  of  the  irregularity  layer  and  the  variation  of  C^, 
or  <An#2>,  within  the  layer. 

The  data  set  from  which  the  estimates  of  these  parameters  are  to 
be  obtained  will  be  collected  by  three  instruments  in  the  DMSP  SSIES 
(Special  Sensor  for  Ions,  Electrons,  and  Scintillation)  sensor 
package.  This  data  set  will  contain  the  following  in  situ 
observations: 

(1)  High  time-resolution  (24  samples/sec)  measurements  of  the  ion 
density  and  measurements  of  the  ion  density  irregularity  PDS  at  high 
fluctuating  frequencies  from  the  ion  Scintillation  Meter  (SM).l3) 

(2)  Measurements  of  the  horizontal  and  vertical  cross-track  ion 
drift  velocities  from  the  ion  Drift  Meter  (DM).(3> 

(3)  Measurements  of  the  ion  and  electron  temperatures,  the 
densities  of  0+  and  the  dominant  light  ion  (H+  or  He+ ) ,  and  the 
horizontal  ram  ion  drift  velocity  from  the  ion  Retarding  Potential 
Analyzer  (RPA).I4! 


*  The  cited  reference  develops  the  theory  in  terms  of  an  earlier 
definition  of  the  strength  parameter,  CG,  which  is  defined  at  a  scale 
size  of  2jt  meters.  It  is  related  to  via  Cs  =  (  in/ 1  000  )3  +  2 


The  basic  data  of  this  set  is  the  high  time-resolution  density  data 
from  the  SM  which  will  be  used  to  generate  estimates  of  the 
irregularity  PDS.  The  drift  velocity  measurements  from  the  DM  and  RPA 
will  be  used  in  calculating  an  estimate  of  C]<  from  parameters  obtained 
from  the  PDS,  and  the  other  measurements  from  the  RPA  will  be  used  in 
calculating  C^L  from  C^. 

The  project  is  divided  into  two  phases.  In  the  first  phase, 
techniques  for  calculating  estimates  of  C^L  from  the  SSIES  data  set 
will  be  developed,  and  parametric  studies  will  be  conducted  to 
determine  the  uncertainties  in  the  final  C^L  estimates  due  to 
uncertainties  in  the  parameters  and  procedures  used  to  calculate  the 
estimates.  There  will  be  no  DMSP  SSIES  data  available  during  this 
phase,  as  the  first  is  due  to  be  launched  in  mid-1987,  so  these 
studies  will  be  made  using  other  data  sources. 

The  second  phase,  which  will  begin  after  the  scheduled  DMSP 
launch,  will  focus  on  how  well  these  techniques  work.  There  will  be 
two  investigations  conducted  during  this  phase:  ( 1 )  an  analysis  of  C^L 
values  calculated  for  selected  DMSP  orbits,  and  (2)  an  assessment  of 
the  validity  of  the  basic  assumptions  made  in  order  to  calculate  an 
estimate  of  C^L  from  a  measurement  made  at  an  altitude  of  830km. 
The  cornerstone  of  the  second  investigation  is  planned  to  be  at  least 
one  coincident  measurement  campaign  during  which  ionospheric  profile 
data  from  an  incoherent  radar,  phase  scintillation  data  from  a 
satellite  beacon,  and  in  situ  ion-density  irregularity  observations 
near  the  F2  peak  would  be  collected  in  near-coincidence  (time  and 
location)  with  a  DMSP  orbit. 


3.  PARAMETRIC  STUDIES:  UNCERTAINTIES  IN  Ck 

The  objective  of  the  first  set  of  parametric  studies  was  to 
determine  the  level  of  uncertainty  in  estimates  of  Ck  from  several 
sources.  Three  studies  were  conducted:  (1)  an  investigation  of 
uncertainties  in  Ck  calculated  from  and  q  due  to  uncertainties  in 
T1  and  q,  (2)  an  investigation  of  uncertainties  in  Ck  calculated  from 
<ANe2>  due  to  uncertainties  in  the  cutoff  frequency,  fc,  and  q,  and 
(3)  an  investigation  of  uncertainties  in  Ck  due  to  uncertainties  in 
the  effective  satellite  velocity,  Vp. 

3.1  Calculation  of  Ck 

According  to  phase-screen  theory,  1*1  an  estimate  for  can  be 
calculated  from  an  In  situ  measurement  of  the  ionospheric  electron 
.density  from 
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Ck  a  5.0x10°  q  ( - )  T1  [1] 

\  VP/ 

where  q  and  T1  are  the  slope  and  intercept  of  a  log-linear  fit  to  the 
Power  Density  Gpectrum  (PDS)  of  the  data  sample,  and  vp  is  the 
effective  velocity  of  the  satellite  with  respect  to  the 
irregularities.  This  last  parameter  is  defined  in  a  coordinate  system 
defined  by  the  irregularities  and  is  given  by 

vp  “  XT  C  v  [2] 


where  v  is  the  vector  velocity  of  the  satellite  and  C  is  a 
transformation  matrix  derived  from  a  generalized  irregularity  model.'5' 
Equation  [  1]  can  also  be  written  in  terms  of  the  variance  of  the 
electron  density  irregularities,  <ANe2>,  as 

q-1 

/  '  ~c  \ 

Ck  =  2.5x10®  q ( q- 1  ) 


<ANe2  > 


[3] 


P 


,v.v. 


where  £c  is  a  low-frequency  cutoff  set  by  whatever  detrending  process 
was  used  in  calculating  <ANe^>,  which  establishes  the  largest  scale 
sizes  included  in  the  variance.  In  this  context,  "detrending" 
includes  all  processing  done  to  the  data  prior  to  calculation  of  the 
variance,  including  selection  of  a  finite-length  data  set. 


3.2  Uncertainties  in  Calculated  from  and  q 

The  effects  of  various  processing  techniques  which  are  used  in 
calculating  the  PD9  were  investigated  using  both  simulated  data  sets 
with  known  spectral  characteristics  and  actual  observations  of 
ionospheric  irregularities.  The  purpose  of  this  study  was  to 
determine  the  optimum  analysis  methodology  for  obtaining  estimates  of 
the  power-law  parameters  q  and  T1  for  use  in  calculating  C^.  This 
study  is  similar  to  a  recent  study  of  the  relative  merits  of  deriving 
estimates  of  the  PDS  from  FFTs  and  from  the  Maximum  Entropy  Method 
( MEM ) , M  but  this  study  will  focus  solely  on  the  FFT  method  and  will 
include  the  effects  of  the  various  processing  methods  on  the  strength 
parameter,  Ti ,  as  well  as  on  the  slope  parameter,  q. 

There  is  currently  no  set  of  data  from  a  DMSP  SSIES  sensor  as  the 
first  of  these  is  scheduled  to  be  on  DMSP  satellite  F8,  due  for  launch 
in  mid-1987.  Two  types  of  data  sets  were  used  in  these  parametric 
studies:  (1)  simulated  density-data  sets  constructed  with  known  power 
density  spectra,  and  (2)  a  sample  of  phase  scintillation  data  taken 
from  the  DNA  WIDEBAND  satellite  experiment.  Details  on  the 
construction  of  the  simulated  data  sets  and  on  the  WIDEBAND  phase 
scintillation  data  are  given  in  Appendix  A  to  this  report. 

All  realizations  for  each  of  the  data  sets,  simulated  and  real, 
were  processed  in  the  following  manner: 

1  .  The  An  data  from  the  simulation  data  sets  were  converted  to 
simulated  plasma  density  samples  using  a  T1  of  1016  and  a  mean  density 
value  of  105  el/cm^.  Optionally,  the  low-frequency  trend  terms  are 
added  on. 


2.  The  plasma  density  sample  (512  points)  is  then  detrended 
using  one  of  three  detrenders  (linear  trend  removal  ( LDET ) ,  quadratic 
trend  removal  ( QDET ) ,  or  end-matching  a  removal  of  the  residual  mean 
( EMMR  )  ) . 

3.  The  detrended  data  are  then  windowed  using  one  of  eleven 
windows  (a  rectangular  window  (essentially  no  window)  or  one  of  ten 
split-bell  cosine  windows  ranging  from  a  10%  taper  to  a  100%  taper). 

4.  An  estimate  of  the  PDS  of  the  sample  is  obtained  from  an  FFT 
of  the  windowed,  detrended  data. 

5.  The  PDS  is  optionally  smoothed  using  one  of  four  moving, 
centered  smoothing  functions  using  binomial  weights  (3-,  5-,  7-,  or  9- 
point  smoothing). 

6.  Estimates  of  T«  and  q  are  obtained  from  a  log-linear  fit  to 
the  PDS  over  the  frequency  range  0.2  to  7.0  Hz. 

The  results  of  this  processing  for  each  realization  ( ,  q,  and 
<ANe^>1/2  at  various  steps  in  the  processing)  are  stored  in  an 
analysis  data  base  for  further  reduction.  A  total  of  1 32  analyses  (3 
detrenders  x  11  windows  x  4  smoothers)  are  stored  in  the  analysis  data 
base  for  each  realization.  Figure  1  shows  an  example  of  the  results 
from  an  analysis  of  the  maximum- leakage  simulation  using  the  quadratic 
detrender,  a  30%  cosine  window,  and  a  3-point  smoother.  The  upper 
plot  shows  the  difference  between  the  measured  slope  (q)  and  the  q 
value  used  to  generate  the  data  set  (q0),  and  the  lower  plot  is  the 
log-difference  between  the  measured  strength  (Tj)  and  the  T,  value 
used  in  the  generation  ( 1 0 1 6  1 .  Both  are  plotted  as  a  function  of  q0. 

The  effects  of  windowing  and  smoothing  on  the  analysis  of  the 
maximum-leakage  simulation  data  set  can  be  seen  in  Figure  2.  These 
plots  show  the  variation  of  Aq  (q  -  q0 )  and  Alog  (^  )  -  16.0)  as  a 
function  oi  percent  window  taper  for  the  non-smoothed  case  and  for 
each  of  the  four  smoothers.  The  low-frequency  terms  were  added  to  all 
data  sets,  and  all  were  detrended  using  a  quadratic  detrender.  Each 
data  point  is  the  average  from  all  66  realizations  for  that  parameter 
and  processing  case.  The  most  striking  effect  shown  in  these  plots  is 
the  systematic  loss  in  spectral  power  (as  measured  by  T-|  )  caused  by 
windowing  the  data  and  the  recovery  of  this  loss  via  the  smoothing 
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figure  2.  Effects  of  windowing  and  smoothing  on  data  from  the  maximum-leakage 
simulation.  All  data  sets  were  detrended  with  a  quadratic  detrender. 
Numbers  Indicate  the  number  of  points  In  the  smoother. 


process.  It  was  expected  that  the  effects  of  windowing  and  spectral 
smoothing  would  be  inter-related,  as  windowing  in  the  time  domain  is 
equivalent  to  smoothing  in  the  frequency  domain,  but  the  loss  of 
accuracy  in  T^  due  to  windowing  and  the  gain  due  to  smoothing  was 
completely  unexpected.  There  is  also  a  gain  in  terms  of  reduction  of 
variance  due  to  the  smoothing.  Figures  3  and  4  show  the  variation  of 
Aq  and  Alog(T^ )  for  the  unsmoothed  and  the  smoothed  (3-point 
smoother)  cases.  The  vertical  bars  indicate  the  standard  deviation 
(a)  of  the  66  analyses  in  each  case.  The  variance  (a2)  within  all 
cases  has  been  reduced  by  roughly  a  factor  of  four  for  both  Aq  and 
Alog  ( T i ) .  (Note:  This  analysis  was  repeated  using  the  linear  and 
end-match/mean-remove  detrenders  with  the  same  results. ) 

Aside  from  this  result,  the  general  effect  of  increased  severity 
in  the  window  was  to  decrease  (though  not  necessarily  improve)  both 
Aq  and  AlogfT^),  and  the  general  effect  of  smoothing  was  to  increase 
(though  not  necessarily  worsen)  both.  The  selection  of  an  optimal 
choice  of  window  severity  and  smoother  size  came  down  to  a  trade-off 
between  accuracy  of  reproducing  q  or  T^ .  The  final  selection  was  to 
use  a  30%  cosine  window  and  a  5-point  smoother.  Figure  5  is  a  similar 
plot  to  Figure  4  for  the  5-point  smoother  showing  the  standard 
deviations  for  this  case. 

Figure  6  indicates  the  effects  of  the  three  detrenders  used  in 
these  tests.  In  both  plots,  the  upper  curve  (labeled  L)  is  from  the 
linear  detrend,  the  next  down  is  from  the  end-match/mean-remove 
detrend  (E),  and  the  lower  curve  from  the  quadratic  detrend  (Q).  As 
can  be  seen,  there  is  little  difference  between  the  three, 
particularly  for  the  more  severe  windows.  For  the  30%  window,  the 
quadratic  detrend  case  is  slightly  better  that  the  other  two,  and  the 
variance  1b  also  slightly  better  for  the  quadratic  detrend. 

The  main  difficulty  in  performing  a  similar  analysis  of  the 
results  of  processing  the  Wideband  data  sample  is  in  determining  what 
"truth"  is  in  terms  of  q  and  T|  for  each  5 1 2-point  data  set.  For  the 
purposes  of  this  study,  it  was  decided  to  define  truth  as  the  results 
obtained  from  the  processing  procedure  selected  in  the  simulation 
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Figure  6.  Detrender  effects  from  analysis  of  the  maximum-leakage  simulation 
data.  A  5-polnt  smoother  was  used  In  all  cases. 
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study  (quadratic  detrend,  30%  cosine  window,  5-point  smoother)  and  use 
these  values  for  calculating  Aq  and  AloglT,).  This  will  at  least 
show  whether  the  general  effects  found  in  the  simulation  study  are 
also  true  for  the  observed  data  set. 

Figure  7  shows  plots  of  Aq  and  Alog  (Tp  as  functions  of  window 
severity  for  each  of  the  five  smoothing  cases  for  pass  PF-52-47.  The 
same  behavior  is  found  in  these  plots  as  in  Figure  2  for  the  simulated 
data  set,  i.e.,  the  loss/gain  in  AlogfTj)  due  to  windowing/smoothing, 
and  the  trends  in  Aq  and  AlogCT^  )  as  functions  of  window  severity  and 
smoother  size.  Note  that  both  curves  pass  through  zero  for  the  30% 
cosine  window  and  5-point  smoother  as  this  was  defined  as  "truth." 
The  variances  for  each  data  point  for  the  various  cases  were  sIbo 
similar  to  those  from  the  simulated  data  set,  although  they  were  more 
strongly  a  function  of  window  and  smoothing,  since  the  base  values  of 
q  and  logC^  )  ware  defined  in  terms  of  the  output  from  one  of  the 
processing  methods. 

The  results  of  the  detrender  study  for  this  data  set  are 
presented  in  Figure  8.  The  fourth  curve  on  these  plots,  labeled  F,  is 
the  results  obtained  using  the  low-pass  filter  detrender  used  in 
generating  the  phase  plots  in  Figure  A-4.  Since  the  quadratic 
detrender  was  used  in  defining  the  base  analysis  values,  the  results 
from  these  plots  do  not  necessarily  indicate  that  the  quadratic 
detrend  results  (curve  Q )  are  better  than  those  from  the  low-pass 
filter  detrend,  but  they  do  indicate  that  the  two  methods  provide 
essentially  the  same  results. 

Based  on  the  results  of  processing  both  the  simulated  and 
observed  data  sets,  the  processing  method  to  be  used  in  calculating 
estimates  of  q  and  T1  from  the  DMSP  SSIES  data  sets  is  as  follows: 

1.  Detrend  each  512-point  data  set.  by  removing  the  quadratic 
trend  determined  by  a  least- squares  fit  to  the  data  set. 

2.  Window  the  detrended  data  using  a  30%  split-bell  cosine 
window . 

3.  Calculate  an  estimate  of  the  PDS  from  an  FFT  of  the  windowt?d, 
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Figure  7. 


Effects  of  windows  and  smoothing  on  data  from  pass  PF-52-47.  All 
data  sets  were  detrended  using  a  quadratic  detrender.  Numbers 
Indicate  number  of  points  In  smont'.er. 
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Figure  a.  Detrender  effects  from  analysis  of  data  from  pass  PF-52-47. 
5-polnt  smoother  was  used  In  all  cases. 


4.  Smooth  the  PDS  estimate  using  a  5-point,  centered  smoother 
with  binomial  weights. 

5.  Calculate  estimates  of  q  and  T,  from  a  log-linear  fit  to  the 
smoothed  PDS  over  the  frequency  range  0.2  to  7.0  Hz. 

The  above  processing  was  performed  on  the  two  simulated  data 
sets,  each  containing  66  data  samples,  and  on  a  third  data  set 
containing  60  data  samples  equally  divided  among  q  values  of  i.o,  2.0, 
and  3.0.  The  RMS  error  in  q  for  the  192  samples  was  3.3%  and  the 
error  in  T1  was  6.7%.  Table  1  shows  the  variation  of  Ack  (%)  with  q 
and  Vp  calculated  using  Equation  [1]  with  these  percent  RMS  errors  in 
q  and  Tj.  The  ACk  values  in  the  bottom  row  are  weighted  RMS  values 
for  each  vp  calculated  using  the  weights  listed  in  the  rightmost 
column.  These  weights  are  generated  from  a  gaussian  function  centered 
at.  q  -  1.8  with  a  1/e-width  of  0.6.  This  is  done  to  simulate  the 
expected  distribution  of  q,  and  was  taken  from  analyses  of  phase 
scintillation  data  collected  from  the  DNA  HiLat  satellite.  The 
expected  errors  in  Ck  due  to  errors  in  the  estimates  of  q  and  T1  range 
from  roughly  7%  at  low  values  of  vp  to  roughly  10%  for  high  values. 

3.3  Uncertainties  in  Ck  Calculated  from  <ANe2> 

The  alternative  method  for  calculating  an  estimate  of  Ck  from  a 
plasma  density  sample,  employed  in  cases  where  the  PDS  is  not 
calculated  from  the  density  samples,  uses  the  density  variance, 
<AN2>.  The  estimate  of  Ck  is  calculated  using  an  Equation  [3]  in 
which  q,  the  PDS  power-law  slope,  is  not  observed  but  rather  set  to 
some  median  value.  As  the  shape  of  the  underlying  PDS  for  ionospheric 
density  irregularities  is  a  red  power-law,  the  variance  for  a  given 
data  sample  1b  very  strongly  dependent  on  both  the  slope  of  the  PDS 
and  the  effects  of  detrending  reflected  in  the  cutoff  frequency,  f  . 
T.he  purpose  of  the  parametric  study  for  this  method  is  to  determine 
the  magnitude  of  Inaccuracies  introduced  in  Ck  calculated  using 
Equation  [3]  due  to  uncertainties  in  the  values  used  for  f  and  q. 


TABLE  1.  Variation  of  Ac^  (*)  with  q  and  v  for  Aq  >  3-3*  and 
ATj  -  6.7*.  - -  **--  ■-  -■  '  *  -  A- 
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1000. 
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4000. 

5000 . 

6000. 
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0.04 
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0.04 
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0.4994 
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0.05 
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6.0 
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1.5 

0.05 
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6.7 
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7.5 
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6.1 

5.5 
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7.4 
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0.8948 

1.7 
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0.9726 

1.8 

0.06 

6.1 
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6.1 
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1.0000 

1.9 

0.06 

6.1 
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7.0 

7.8 

8.5 
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0-9726 

2.0 

0-07 

61 
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0.07 
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6.5 
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0.7788 
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0.07 
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0.08 
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0.3679 

6.1 

5.5 

6.1 
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Estimates  for  fc  and  the  uncertainty  in  this  parameter  for  the 
quadratic,  linear,  and  end-match  detrenders  were  obtained  by 
calculating  an  effective  fc,  defined  by 

2T1 

Ce  (q-1 )<AN2> 

for  each  simulated  realization  and  calculating  the  mean  and  variance 
of  the  values  obtained  for  each  detrender.  This  provided  the 
following  results: 

Quadratic  detrend:  fc  -  0.06  +-25% 

Linear  detrend:  fc  ■  0.04  +-45% 

End-match/mean  remove:  fc  ■  0.03  +-50%. 

The  value  of  fc  for  the  low-pass  filter  detrender  is  essentially 
determined  by  the  cutoff  frequency  specified  in  constructing  the 
filter.  For  the  sake  of  this  study,  we  will  define  the  detrender 
cutoff  frequency  such  that  the  detrend  period,  which  is  the  reciprocal 
of  this  frequency,  is  roughly  one-half  the  sample  interval.  Since  the 
simulation  data  are  set  up  to  be  24  samples/second  and  the  sample  size 
is  512  data  points,  the  sample  interval  is  21.333...  seconds.  For 
convenience,  we  will  select  a  10-second  detrend  interval  which 
provides  fc  ■»  0.10.  It  was  difficult  to  determine  a  variance  in  fc 
for  this  detrender  using  the  simulated  data  sets  because  these  data 
were  constructed  in  such  a  way  that  the  power  is  located  at  discrete 
frequencies  rather  than  spread  across  all  frequencies.  In  the  tests 
run,  however,  the  largest  variances  found  were  on  the  order  of  3-4%. 

An  initially  attractive  feature  of  this  method  for  calculating  Ck 
is  that  one  need  not  calculate  the  PDS  of  the  data  sample  if  a  mean 
value  of  the  PDS  slope,  q,  could  be  used.  In  early  simulations, 
however,  it  became  apparent  that  this  could  not  be  done.  Even  for  the 
best  cases,  assuming  no  error  in  fc  or  any  other  parameters,  the 
minimum  RMS  errors  in  were  on  the  order  of  50%.  The  RMS  error 
increased  to  over  100%  when  realistic  values  for  expected  errors  In 


other  parameters  were  used.  It  was  decided,  therefore,  to  assess  the 
level  of  errors  to  be  expected  from  this  method  assuming  that  q  is 
measured  with  a  level  of  accuracy  established  in  earlier  studies  (+- 
3.3%). 

Tables  2a-c  summarize  the  expected  errors  in  using  Equation  2 
for  the  Fourier  detrender  ( fc  «  0.10,  Afc  =  3.5%),  the  quadratic 
detrender  ( f c  ■  0.06,  Afc  >  25%),  and  the  linear  and  end-match 
detrenders  ( fc  ■  0.35,  Afc  ■  45%).  The  errors  for  the  Fourier 
detrender  are  comparable  to  those  found  for  the  calculation  method 
using  q  and  T1  shown  in  Table  1 ,  while  the  levels  for  the  quadratic 
and  linear  detrenders  are  factors  of  2  and  4  larger,  respectively. 
[Note:  As  with  the  results  reported  earlier  for  q  and  T1 ,  these  tables 
do  not  include,  the  effects  of  expected  errors  in  vp.]  The  large 
errors  in  the  quadratic  and  linear  detrender  cases,  reflected  in  the 
large  variances  in  the  effective  values  for  fc,  are  due  to  the  data- 
dependent  nature  of  the  effects  of  these  detrenders  on  the  large-scale 
features  which  dominate  <AN2>.  In  other  words,  the  amount  of  low- 
frequency  power  removed  in  the  detrending  process  will  depend  on  the 
location  of  the  extrema  of  the  large-scale  features  in  the  data 
sample.  The  Fourier  detrender,  which  must  be  run  on  the  entire  data 
set  prior  to  selecting  512-point  data  samples  for  processing,  is 
designed  to  remove  power  only  at  frequencies  below  the  specified 
cutoff  frequency  and  is  not  affected  by  the  location  of  the  extrema. 

In  summary,  if  this  method  is  to  be  used  for  calculating  C^,  the 
following  rules  should  be  followed: 

a.  The  data  should  be  detrended  using  a  detrender  which  will 
remove  power  at  the  larger  scales  in  a  predictable  manner,  such  as  the 
Fourier  detrender  used  in  this  study. 

b.  A  power-density  spectrum  ( PDS )  should  be  constructed  from  each 
data  sample  as  described  in  Section  4.1  to  obtain  an  estimate  of  the 
spectral  slope  rather  than  using  a  mean  value  for  q. 

If  the  data  are  processed  in  this  fashion,  the  expected  errors  in  0^ 
for  this  technique  due  to  uncertainties  in  fc  and  q  will  be  roughly  6- 
9%  for  low  values  of  vp  and  12%  for  high  values. 
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TABLE  2c-  Variation  of  Ac^  (H)  with  q  and  vp  for  the  linear 

detrenders  (fc«0035t  Afc-45*)  using  Aq  >  3.3*.  The 
bottom  line  Is  the  weighted  RMS  Ack  for  ell  q  as  a 
function  of  Vp. 
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3.4  Uncertainties  in  Due  to  Uncertainties  in  Vp 

The  final  part  of  this  parametric  study  is  aimed  at  determining 
the  expected  errors  in  estimating  vp  and  the  effects  of  these  errors 
on  estimates  of  C^.  This  parameter  is  a  function  of  (1)  the  satellite 
velocity  vector,  which  is  well  determined;  (2)  the  orientation  of  the 
geomagnetic  field  at  the  satellite,  which  can  be  fairly  well 
determined  from  standard  models  of  the  geomagnetic  field,  (3)  the 
velocity  of  the  irregularities  with  respect  to  the  satellite,  which 
can  be  measured  or  taken  from  a  model;  and  (4)  the  axial  ratios  of 
the  irregularities  and  their  orientation  with  respect  to  the 
geomagnetic  field,  which  cannot  be  extracted  from  the  satellite  data 
set  and  must  be  obtained  from  models  which  are  still  in  their  early 
evolutionary  stages.  This  study  will  focus  on  the  effects  of  errors 
in  the  measured/modeled  drift  velocity  and  in  the  modeled  axial  ratios 
on  Vp,  and  the  effects  of  those  errors  on  Ck. 

Figure  9  is  a  scatter  plot  of  vp  vs.  apex  (magnetic)  latitude 
calculated  for  28  simulated  DMSP  orbits,  14  of  them  with  the  local 
time  of  ascending  node  at  0600  and  14  at  1000.  The  IGRF80  model  was 
used  to  obtain  parameters  pertaining  to  the  the  geomagnetic  field,  and 
modified  WBMOD  models  were  used  for  the  axial  ratios  (a  and  b)  and  the 
in  situ  drift  velocity.  The  pattern  in  this  plot  is  a  function  of  the 
orientation  of  the  DMSP  orbit  with  respect  to  the  geomagnetic  field 
and  the  orientation  and  axial  ratios  of  the  irregularities.  In  the 
equatorial  region,  the  velocities  are  low  because  (i)  the  angle 
between  the  orbital  plane  and  the  geomagnetic  field  direction  is  small 
and  (2)  the  irregularities  in  the  equatorial  region  are  elongated  rods 
aligned  with  the  geomagnetic  field  direction.  Thus  the  satellite  is 
moving  along  the  long  axis  of  the  irregularities,  resulting  in  a  low 
Vp.  The  increase  in  vp  with  apex  latitude  is  due  primarily  to  the 
changing  angle  between  the  satellite  velocity  vector  and  the 
geomagnetic  field  vector.  The  sudden  "spreading"  of  the  vp  curve  at 
high  latitudes  is  due  to  (  I  )  the  addition  of  an  appreciable  in  situ 
drift  velocity  as  the  satellite  enters  the  high- latitude  convection 
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pattern  and  (2)  the  introduction  of  sheet-like  irregularities  (b>i)  in 
certain  sections  of  the  auroral  zone.  [Note:  In  the  model  used  for 
axial  ratio  b,  sheet-like  irregularities  are  confined  to  the  evening 
side  of  the  auroral  oval  (roughly  1800-2400  magnetic  local  time 
(MLT)).]  As  we  shall  see  later,  the  effect  of  the  sheet-like 
irregularities  is  to  reduce  the  calculated  vp. 

In  order  to  estimate  the  effects  of  errors  in  a  and  b,  vp  was 
calculated  for  the  dawn-dusk  (0600LT)  cases  used  in  generating  Figure 
9  with  induced  "errors"  of  +-50%  in  each  axial  ratio,  and  the  percent 
error  in  vp  was  calculated.  The  results  are  shown  in  the  two  scatter 
plots  in  Figure  10.  As  can  be  seen  in  the  upper  plot  in  this  figure, 
Vp  is  very  insensitive  to  errors  in  a  except  for  a  few  cases  near  the 
magnetic  equator.  All  but  3  of  1440  data  points  resulted  in  an  error 
in  vp  of  less  than  5%,  and  the  maximum  error  in  vp  was  26%.  On  the 
other  hand,  vp  can  be  very  sensitive  to  errors  in  axial  ratio  b,  as 
can  be  seen  in  the  lower  plot  in  Figure  10.  [Note:  Errors  were 
introduced  only  in  those  cases  where  b>i,  and  b  was  not  allowed  to  be 
less  than  1.]  The  bulk  of  the  errors  is  still  below  5%  (96%  of  the 
errors  at  all  latitudes,  and  88%  of  the  errors  for  latitudes  >60 
degrees,  were  below  5%);  however,  the  range  in  vp  errors  is  now  -35% 
to  +100%.  Fortunately,  the  area  where  sheet-like  irregularities  exist 
appears  to  bo  limited  to  the  evening-side  auroral  oval  (Fremouw, 
private  correspondence);  unfortunately,  this  is  an  area  of  interest. 

A  second  error  which  involves  the  axial  ratios  is  specifying 
sheet-like  irregularities  where  the  irregularities  are  actually  rod¬ 
like,  and  vice  versa.  Figure  1 1  is  a  scatter  plot  of  vp  percent  error 
for  the  same  cases  used  in  Figure  10  when  b  is  set  to  1.0  everywhere. 
The  range  of  errors  is  now  0%  to  over  300%,  with  95%  of  the  total 
distribution,  and  84%  of  the  distribution  at  latitudes  >60  degrees, 
with  errors  <5%.  Figure  12  shows  the  effect,  of  forcing  rod  like 
irregularities  everywhere  on  vp.  The  upper  plot  shows  vp  on  lml.it  .>d 
using  the  model  for  b  for  the  dawn-dusk  orbits,  and  tin-  lower  plot 
shows  vp  calculated  for  1 .  Removing  the  sheet  like  irregularities 
shifts  all  values  of  vp  above  600U  m/s.  Again,  the  locations  and 


Figure  10.  Effects  of  ±50%  errors  in  axial  ratio  a  (upper  plot)  and 
b  (lower  plot)  on  parameter  Vp. 
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times  where  the  largest  errors  occur  will  be  limited,  but  they  are  in 

latitude  regimes  of  Interest. 

A  second  potential  source  of  error  in  Vp  is  errors  in  the 
measured,  or  modeled,  values  of  the  in  situ  drift  velocity.  Figure  13 
shows  the  effect  on  Vp  of  removing  the  drift  velocity  from  the 
calculation.  The  errors  range  from  +-25%  at  high  latitudes  to  +-100% 
at  equatorial  latitudes.  The  overall  RMS  error  was  10%,  with  16%  at 
equatorial  latitudes  and  6%  at  high  latitudes.  The  larger  errors  at 
equatorial  latitudes  are  due  to  the  east-west  component  of  the  drift 
velocity,  which  changes  the  angle  between  the  vector  velocity  of  the 
irregularities  with  respect  to  the  satellite  and  the  local  geomagnetic 
field  vector.  When  this  angle  is  small  to  begin  with,  as  is  the  case 
for  DM8P  orbits  in  the  equatorial  region,  Vp  is  very  sensitive  to  this 
angle. 

The  effects  of  errors  in  the  horizontal  drift  velocities  An  the 
across-  (u^)  and  along-  (ur)  orbit  track  are  shown  in  Figure  14.  The 
upper  plot  shows  the  errors  introduced  by  +-50%  errors  in  ur  and  the 
lower  plot  for  +-50%  errors  in  u^.  [Note:  Errors  in  vp  introduced  by 
errors  in  the  vertical  component  of  the  drift  velocity,  uv,  were 
routinely  <1%.]  The  errors  induced  in  vp  by  errors  in  ur  are  small  at 
equatorial  and  mid  latitudes  (<1%),  and  are  bounded  by  +-15%  at  high 
latitudes.  The  errors  induced  by  errors  in  u^  are  smaller  at  high 
latitudes  ( <  1 0 %  ) ,  but  can  be  much  larger  in  the  equatorial  region  (up 
to  60%)  for  the  reasons  discussed  in  the  previous  paragraph.  The 
overall/equatorial/high  latitude  RMS  values  for  the  two  cases  were 
7%/9%/2%  for  errors  in  u^  and  2%/<l%/4%  for  errors  in  ur. 

Although  it  is  relatively  straightforward  to  estimate  the 
propagation  of  errors  from  the  axial  ratios  or  the  drift  velocities  to 
Vp,  it  is  not  so  simple  to  make  estimates  of  what  errors  are  to  be 
expected  in  the  parameters  themselves.  Both  axial  ratios  must  be 
provided  from  a  model  (taken  from  the  WBM0D  scintillation  model  and 
slightly  modified)  which  is  based  on  values  inferred  from  observations 
of  phase  scintillation  rather  than  direct  observations  of  these 
parameters.  While  it  appears  errors  in  axial  ratio  a  will  not 
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drastically  effect  vp,  errors  in  axial  ratio  b  can  cause  errors  in  vp 
greater  than  100%.  Unfortunately,  there  is  more  confidence  in  the 
model  for  a  than  for  b,  and  no  good  estimate  of  how  "bad"  the  model 
for  b  is.  The  estimate  of  +-50%  used  in  the  parametric  study  is 
probably  as  good  as  any  at  this  time. 

The  outlook  is  not  quite  so  bleak  for  estimating  errors  in  the 
drift  velocities,  assuming  that  drift  velocity  measurements  are 
available  from  the  SSIES  Drift  Meter  (DM)  and  Retarding  Potential 

Analyzer  ( RPA )  sensors.  The  DM  sensor  will  provide  estimates  of 

and  uv,  and  the  RPA  will  provide  ur.  In  parametric  studies  conducted 
while  developing  software  to  process  these  data,t?l  it  was  found  that 
the  probable  minimum  RMS  error  in  these  measurements  was  30  m/s  for 
the  DM  measurements  and  100  m/s  for  the  RPA  measurement.  Figure  15 
shows  the  percent  error  in  vp  for  the  dawn-dusk  orbits  for  errors  of 
+-50  m/s  in  u^  and  uv  and  +-150  m/s  in  ur.  The  +-2%  offsets  shown  in 

the  figure  are  due  to  the  errors  in  ur,  and  the  scatter  at  equatorial 

and  high  latitudes  is  due  to  errors  in  u^.  The  RMS  error  at  all 
latitudes  is  2%,  dominated  by  the  error  in  ur,  and  the  99th  percentile 
is  about  6%. 

If  the  drift  velocities  are  not  available,  however,  errors  on  the 
order  of  those  shown  in  Figure  13  must  be  accepted  or  a  model  for  the 
drift  velocity  must  be  used.  While  it  is  unlikely  that  a  simple  model 
will  provide  accurate  drift  velocity  measurements  at  high  latitudes, 
it  may  be  possible  to  model  the  equatorial  east-west  drifts  to  within 
50%.  Figure  16  shows  the  expected  errors  in  vp  for  the  case  where  no 
drift  velocities  are  used  at  high  latitudes,  and  a  model  with 
accuracies  of  +-50%  is  used  at  equatorial  latitudes.  The  RMS  errors 
in  vp  for  this  case  are  7%  for  all  latitudes,  9%  at  equatorial 
latitudes,  and  8%  at  high  latitudes  with  99th  percentile  values  of 
about  25%. 

In  summary,  the  expected  errors  in  vp  are  as  follows: 

a.  If  reasonably  good  observations  of  the  in  situ  drift 
velocities  are  available,  errors  in  vp  will  be  on  the  order  of  a  few 
percent  at  equatorial  and  mid  latitudes  and  at  high  latitudes  where 
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Figure  15.  Effect  of  expected  minimum  errors  In  drift  velocity  measure 
ments  on  parameter  V  . 


the  irregularities  are  rod-like.  At  high  latitudes,  where  they  may  be 
sheet-like  (most  probably  in  the  evening-side  of  the  auroral  oval), 
the  RMS  errors  will  be  on  the  order  of  15-25%,  with  errors  greater 
than  100%  not  uncommon,  due  strictly  to  uncertainties  in  (1)  where  to 
switch  from  rods  to  sheets  and  (2)  the  value  to  use  for  axial  ratio  b 
for  sheet-like  irregularities. 

b.  If  observations  are  not  available,  the  expected  RMS  errors  in 
both  the  equatorial  region  and  the  high  latitude  region  (away  from  the 
region  of  sheet-like  irregularities)  increase  to  about  10%  due  to 
uncertainties  in  the  drift  velocities.  Errors  in  the  high  latitude 
region  where  sheet-like  irregularities  may  be  found  will  probably 
increase  to  20-30%. 

Figure  17  closes  the  loop  on  this  part  of  the  study, 
demonstrating  the  expected  RMS  error  in  for  a  given  RMS  error  in 
Vp.  The  RMS  AC)<  values  were  calculated  using  the  assumed 

distribution  of  q  used  in  the  previous  calculations  (a  gaussian 
distribution  centered  at  1.8,  with  1/e-width  of  0.6).  In  the  range  - 
20%  to  +20%,  the  relationship  between  the  two  is  nearlv  linear,  so  the 
expected  errors  in  vp  described  in  the  previous  paragraph  are  fairly 
good  estimates  for  the  expected  RMS  errors  in  due  to  orrora  in  vp. 

3.5  Summary  of  Cj^  Uncertainties 

In  the  preceding  sections  we  found  that  the  expected  errors  in 
due  to  processing-induced  errors  in  q  and  T1  or  in  q  and  fc  are  on  the 
order  of  5-10%,  and  the  errors  due  to  errors/uncertainties  in  the 
calculated  value  of  vp  range  from  a  few  percent  to  over  100%.  It  is 
not  surprising,  therefore,  to  find  that  the  controlling  factor  on  how 
accurately  can  be  measured  is  the  uncertainties  in  vp.  This  can  be 
seen  in  Tables  3a  and  3b,  which  summarize  the  expected  uncertainties 
in  as  a  function  of  vp  and  Avp  for  each  of  the  two  methods  of 

calculating  C^.  When  the  uncertainties  in  vp  are  low,  J.e.  when  the 
axial  ratios  and  in  situ  drift  velocities  are  determined  accurately, 
the  uncertainties  in  will  be  on  the  order  of  5-15%.  This  will 
usually  be  the  case  In  the  equatorial  regions.  As  the  uncertainties 
in  Vp  increase,  however,  the  uncertainties  in  will  increase  roughly 
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Figure  17.  Effect  of  errors  In  parameter  Vp  on  C^. 
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linearly  in  Avp.  This  will  occur  in  those  regions  where  the  shape  of 
the  irregularities  is  in  question,  i.e.  whether  they  are  rod-like  or 
sheet-like,  and  where  the  in  situ  drift  velocities  are  not  well 
observed  or  modeled.  These  uncertainties  will  be  greatest  at  high 
latitudes,  particularly  (1)  during  geomagnetic  disturbances,  when  the 
drift  velocities  become  difficult  to  measure  or  model,  (2)  in  and 
near  the  evening  sector  of  the  auroral  oval,  where  sheet-like 
irregularities  are  found,  and  (3)  near  the  dayside  cusp  and  nightside 
exit  regions  of  the  high-latitude  convection  pattern  and  the  velocity 
shear  regions  near  the  boundaries  of  the  high-latitude  current 
systems,  where  the  drift  velocities  can  change  dramatically  both  in 
space  and  time. 
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4.  PARAMETRIC  STUDIES:  UNCERTAINTIES  IN  C^L 

The  objective  of  the  second  phase  of  parametric  studies  was  to 
determine  the  level  of  uncertainty  in  estimates  of  CkL  made  from  the 
estimate  of  Ck  at  an  altitude  of  (nominally)  830km.  Two  studies  were 
planned,  one  focusing  on  the  effects  of  uncertainties  in  the  variation 
of  the  background  ionospheric  electron  density  with  height,  and  the 
second  on  the  effects  of  uncertainties  in  the  variation  of  the 
irregularity  distribution  with  altitude.  The  first  study  was 
completed  during  the  period  covered  by  this  report,  the  second  will  be 
completed  during  the  next  contract  year. 


4.1  Calculation  of  C^L  From  Ck 

In  order  to  calculate  an  estimate  of  CkL  from  an  observation  of 
Ck  at  some  altitude,  h8,  the  following  assumptions  must  be  made: 

( 1 )  The  geometry  of  the  irregularities  and  Lhe  slope  of  the  PDS 
of  the  irregularities  are  relatively  constant  throughout  the 
irregularity  layer. 

(2)  The  variation  of  Ck,  or  <ANe2>,  with  altitude  is  known. 

(3)  The  altitude  at  which  the  measurement  is  made  is  within  the 
irregularity  layer. 

(4)  The  measured  Ck  is  representative  of  conditions  throughout 
the  irregularity  layer. 

(For  the  time  being,  we  will  beg  the  issue  as  to  whether  any  of  these 
assumptions  are  valid  or  warranted  but  will  accept  them  as  being 
reasonable.)  Assuming  that  the  irregularity  PDS  and  geometry  remain 
fairly  constant  throughout  the  irregularity  layer,  Ck  can  be  related 
to  CkL  through 


CkL  = 


where  and  are  the  top  and  bottom  of  the  irregularity  layer. 
Since  Ck  is  proportional  to  <ANe2>,  this  can  be  rewritten 
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where  C^thj)  and  <ANe2>8  are  measured  at  the  satellite  altitude.  At 
this  point  some  model  for  the  height  variation  of  <ANe2>  is  roquired. 

The  model  for  <ANe2>  must  specify  (1)  the  height  range  over 
which  irregularities  are  found  (h^  and  h^. ) ,  and  (2)  the  variation  of 
<An82>  over  that  range.  A  sample  model  would  be  to  assume  that  the 
irregularity  layer  extends  from  hmF2  to  the  satellite  height,  hs,  and 
that  <AN82>1/2/N8  remains  constant  throughout  the  layer.  For  this 
case,  Equation  [6]  becomes 


Ne2  ( h )  dh 

V2 _ 

Ne2(hs) 


Ck(h8) 


which  depends  solely  on  the  topside  electron  density  profile  between 
hmF2  and  the  satellite  altitude.  While  this  sole  dependence  on  NtJ(h) 
is  due  to  the  assumed  variation  of  <ANe2>  with  height,  it  is  not 
unreasonable  to  assume  that  the  conversion  of  to  C^L  will  depend  in 
some  way  on  the  topside  electron  density  profile. 

An  effective  irregularity-layer  thickness  parameter,  Le£f,  can  be 
defined  as  the  ratio  of  C^L  to  C)<(h6).  From  Equation  [7],  valid  tor 
the  case  <ANe2>'/2/Ne  ■  constant,  is  then 


I 

i 


t  n< 

nmF2 


2(h)  dh 


Ne2(hs) 


The  calculation  of  C^L  from  is  then  just 


CrL  ■  Cjc(hg)Laj£ 


[9] 


For  a  simple  Chapman-layer  topside  with  a  constant  scale  height,  the 
integral  in  Equation  [8]  can  be  evaluated  analytically  to  provide 


Leff  •  >H 


Ne(hs) 


[10] 


where  -y  is  a  function  of  the  height  of  the  satellite  above  the  F2 
peak  and  the  scale  height  given  by 


us  +  1 


r(ug-2) 


3 

5 


ug  ■  2  exp  < - 


<hg-hrnF2) 


! 


and  NmF2  is  the  density  at  the  F2  peak. 

The  form  of  Equation  [io]  suggests  the  definition  of  a  normalized 
effective  layer-thickness,  in  which  the  ratio  of  the  electron 
density  at  the  satellite  to  that  at  the  peak  is  removed.  This 
parameter  would  then  be  defined  by 


<Leff>  “ 


Ng ( hg ) 

V2 


Jef  f ' 


[n] 


and  the  Equation  [9]  could  be  rewritten 

|2 


Cj^L  ■*  C)<(  hB 


MV 

"m« 


<Lef f > • 


[12] 


[Note  that  for  the  Chapman  profile  case,  <  L  e  f  f  >  is  just  *>H.] 
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4.2  Uncertainties  Due  to  Profile  Model 

The  study  of  the  effects  of  errors  in  parameters  which  define  the 
topside  electron  density  profile  will  focus  on  the  variation  of  the 
two  layer-thickness  parameters,  as  defined  by  Equations  [7]  (LQ£j)  and 
[11]  (<Leff>),  as  functions  of  f0F2,  hmF2,  Yt ,  hT,  and  Na(hs).  The 
profile  model  described  in  the  previous  section  of  this  report,  and 
the  <ANa2>1/2/N#  ■  constant  model  for  the  height  variation  of  <ANe2>, 
were  used  for  this  study. 

Rather  than  look  at  the  variation  of  errors  all  along  a  DMSP 
orbit,  it  was  decided  to  focus  on  those  latitudes  and  local  times 
along  the  orbits  at  which  measurable  levels  of  scintillation  can  be 
expected.  Table  4  lists  the  latitudes  and  local  times  selected 
(basically  equatorial  at  2200LT,  and  auroral  at  all  local  times),  the 
f0F2,  hmF2,  Y{.,  and  hT  values  used  to  construct  the  base  profiles,  and 
the  values  for  and  <Lajf>  for  the  base  profiles.  (Note: 
Initially  a  much  larger  number  of  profiles  were  to  be  processed,  but 
the  results  were  not  much  different  from  profile  to  profile  and  are 
well  represented  by  the  profiles  listed  in  Table  4.)  The  variation  of 
L'eff  and  <Lajf>  with  each  parameter  listed  above  was  calculated  as 
follows: 

(1)  A  base  profile  was  constructed  using  f0F2,  hmF2,  Yt,  and  hT 
from  Table  4.  Values  for  Na(hs),  Laff,  and  were  calculated  for 
this  profile. 

(2)  The  parameter  selected  was  varied  from  -20%  to  +20%  of  the 
initial  value  in  5%  steps. 

(3)  The  profile  was  adjusted  to  fit  the  data  set  consisting  of 
f0F2,  hmF2,  Y^.,  h<p,  and  Ne(h8)  by  iteratively  changing  the  Mg  and  li 
profile  parameters  until  Na(hs)  calculated  from  the  new  profile  was 
within  0.1%  of  the  desired  value. 

(4)  Values  for  LQj£  and  <  Le  f  f  >  were  calculated  for  the  new 
profile,  and  the  percent  change  from  the  corresponding  values  for  the 
base  profile  were  calculated. 


& 


TABLE  <t  •  Profile  parameters  for  L#£f  parameteric  studies  and 

the  values  for  Lg££  and  <L#££>  for  each  basic  profile. 


Each  of  the  24  profile  sets  listed  in  Table  4  were  processed  in  this 
fashion,  and  statistics  for  the  percent  variation  in  Lejf  and  <Leff> 
were  generated  for  the  Auroral  and  Equatorial  cases  separately  and 
together.  Figure  18  shows  samples  of  the  basic  profile  (solid  line) 
and  the  +20%  profiles  (dotted  lines)  for  all  parameters.  The  profile 
case  used  is  labeled  in  each  plot. 

Summaries  of  the  results  are  given  in  Table  5  for  the  Auroral 
cases,  Table  6  for  the  Equatorial  cases,  and  Table  7  for  all  cases 
combined.  The  tables  list  the  average,  variance,  maximum,  and 
minimum  values  for  and  A<LSf{>  (%)  for  the  aggregate  set  under 
each  category.  The  results  for  shown  in  these  tables  can  be 
summarized  as  follows: 

(1)  f qF2 .  Uncertainties  in  this  parameter  resulted  in  the 
largest  corresponding  uncertainties  in  L0jf,  with  a  range  of  55%  to 
+  100%.  This  is  unfortunate,  as  there  will  probably  be  no  direct 
observation  of  this  parameter  over  most  of  the  DMSP  orbit. 

(2)  hmF2.  Although  the  changes  to  the  profile  shape  are  the 
most  dramatic  for  this  parameter  (see  Figure  18),  the  effect  of 
uncertainties  in  hmF2  on  Lg!ff  are  much  smaller  than  the  effects  of 
uncertainties  in  either  f0F2  or  Ne(hg).  The  maximum  error  in  was 
+17%,  obtained  from  a  (probably)  pathological  case  with  high  hmF2  and 

values  and  a  low  hp  value  (E-19).  With  the  exception  of  this  one 
case,  the  maximum  values  for  ALe££  were  all  less  than  10%. 

(3)  h<p.  The  effect  of  uncertainties  in  hip  on  Le£f  was  largely  a 
function  of  whether  hT  was  below  or  Just  above  the  height  of  the 
satellite.  The  effects  were  very  small  for  all  Auroral  cases  for 
which  hp  was  set  to  1200km,  and  were  large  only  for  those  Equatorial 
cases  for  which  hp  was  below  1000km.  Even  for  these  cases,  however, 
the  maximum  error  in  Le££  due  to  uncertainties  in  hT  was  only  +10%  for 
+  20%  errors  in  hp.  On  average,  the  expected  errors  in  Leff  due  to 
uncertainty  in  hp  should  not  exceed  a  few  percent. 

(4)  Yf  This  parameter  was  the  only  one  to  show  a  significant 
and  systematic  difference  between  the  Auroral  and  Equatorial  cases. 
For  the  Auroral  cases,  the  effects  of  uncertainties  in  were 
comparable  to  those  due  to  uncertainties  in  hmF2.  with  a  range  of  -6% 


AUttede  (ka)  _  _  Altitude  (ka) 


k? 
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4N#<hj) 
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log  (N#)  (ol/em5) 


Figure  18.  (Continued). 
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TABLE  5a.  Summary  of  results  for  parameter  LQff  (Auroral). 


Average  &L#ff  (*) 


PARAMETER 

-20* 

-15* 

-10* 

-5* 

♦  5* 

+  10* 

♦  15* 

+  20* 

f_F2 

-52.4 

-42. 0 

-30.0 

-16.0 

18.3 

39-0 

62.6 

89-1 

h°F2 

6.9 

5.2 

3.4 

1.7 

-1.7 

-3.3 

-4.8 

-6.2 

hj 

-2.5 

-1.6 

-0-9 

-0.4 

0.3 

0.5 

0.7 

0-8 

1 

Yt 

-6.3 

-4.9 

-3.3 

-1.7 

1.7 

3.4 

5.4 

7.3 

N,(hs) 

47.1 

32.4 

19.9 

9.2 

-8.1 

-15.0 

-21.) 

-?6. 6 

RMS  Al-££ 

(*) 

PARAMETER 

-20* 

-15* 

-10* 

-5* 

♦  5* 

+  10* 

♦  15* 

X  1 

O  1 

N  1 
+  1 

fnF2 

1.4 

1.2 

0.9 

0-5 

0.6 

1.3 

2.2 

3.2 

h°F2 

1.9 

1.3 

0.7 

0.3 

0.3 

06 

0-9 

1.2 

in 

Ht 

0-3 

0.2 

0.1 

0.1 

0.1 

0.1 

0.1 

0.2 

1 

Yt 

1.2 

0.9 

0.6 

0-3 

0.3 

0.7 

0.9 

1.3 

W 

1.6 

1.1 

0.6 

0.3 

0.2 

0-5 

06 

0.8 

R 


V 


¥; 

I 


V, 

y 

v 

V 


>• 

K 


PARAMETER 


W 


-50.1 

2.4 

-1.9 

-4.7 

45.0 


Minimum  ( *) 


-15.3 

0.9 

-03 

-1.2 

8 


Maximum  Ai,off  (*) 


59.8 

-3.2 

0.4 

4.0 

-20? 


85.0 

-4.1 

0-5 

5-5 

-25-4 


PARAMETER 

-20* 

•  15* 

-10* 

-5* 

♦  5* 

♦  10* 

♦  15* 

♦  20* 

f-F2 

-54.5 

-43-8 

-31.4 

-16.8 

19.2 

40.9 

65-8 

93.7 

h_F2 

9.9 

7.2 

4.6 

2.3 

-2.2 

-4.4 

-6.4 

-8.5 

Ul 

hw 

-3.0 

-1.9 

-1.1 

-0.5 

0.4 

0.7 

0.9 

1.0 

i 

Yt 

-9.0 

-6.8 

-4.8 

-2.4 

2.3 

4.9 

7.5 

10.1 

w 

49.4 

34.0 

20-8 

9.6 

-8.4 

-15-7 

-22.1 

-27.8 

■!; 

•> 

v 

<>■ 

~V', 

$■ 

£ 

£ 

y 

y 

I*1 

y 
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[ABLE  3b.  Summary  of.  results  for  parameter  <laff>  (Auroral). 


f 

! 


Average  A<Lgff>  (*) 


PARAMETER 

-20* 

-15* 

-10* 

f0F2 

16.4 

11.1 

6.8 

V2 

7.0 

5.2 

3.5 

hT 

-2.5 

-1.6 

-0.9 

Tt 

-6  >  3 

-4.8 

-3.3 

N,(Kb) 

-5-7 

-4.3 

-2.9 

RMS 


PARAMETER 

-20* 

-15* 

-10* 

f0F2 

3-4 

2.2 

1.3 

h„F2 

1.9 

1.2 

0.7 

hT 

0-3 

0.2 

0.1 

1.2 

0.9 

0.6 

Nt(h,) 

1.0 

0.8 

0.5 

-3* 

♦  5* 

♦  10* 

♦  13* 

♦  20* 

3.1 

-2.7 

-5.0 

-7.0 

-8.8 

1.7 

-1.6 

-3.2 

-4.7 

-6.2 

-0.4 

0.3 

0.5 

0.7 

0.8 

-1.7 

1.7 

3.5 

5.4 

7.3 

-1.4 

1.5 

2.9 

4.3 

5.8 

A<L,£f> 

-5* 

(*) 

♦  5* 

♦  10* 

♦  15* 

♦  20* 

0.6 

0-3 

0.9 

1.2 

1.3 

0.3 

0.3 

0.6 

0.9 

1.2 

0.1 

0.1 

0-1 

0-1 

0.2 

0.3 

0.3 

0.6 

0.9 

1.3 

0-3 

0.3 

0-5 

0.8 

1.1 

a 


PARAMETER 


f0F2 

y2 

hij« 


w 


PARAMETER 


fcF2 

y2 

nT 

Ng(hg) 


Minimum  (*) 


-20* 

-15* 

-10* 

-3* 

♦  5* 

♦  10* 

♦  13* 

♦  20* 

11. i 

7.6 

4.7 

2.2 

-1.9 

-3.6 

-5.2 

-66 

2.4 

2.2 

1.7 

0.9 

-1.0 

-2.0 

-3-0 

-4.1 

-1.9 

-1.2 

-0.7 

-0.3 

0.2 

0  4 

0.5 

0.5 

-4.7 

-3.6 

-2.5 

-1.3 

1.3 

2.7 

4.0 

5.5 

-4.2 

-3.1 

-2.1 

-1.0 

1.0 

2-0 

3.0 

4.0 

Maximum 

A<Lg££>  (*) 

-20* 

-15* 

-10* 

ar  t 
aft  1 
i  l 
i 

♦  3* 

♦  10* 

♦  13* 

♦  20* 

21.8 

14.6 

8.8 

4.0 

-3.3 

-6.2 

-8.6 

-10.8 

9,9 

7.2 

4.7 

2.3 

-2.2 

-4.3 

-6.4 

-8-5 

-3.0 

-1.9 

-1.1 

-0.5 

0.4 

0-6 

0.9 

1.0 

-9.0 

-6-8 

-4.6 

-2.4 

2.4 

4.9 

7.5 

10. 1 

-7.1 

-5.4 

-3.6 

-1.8 

1.9 

3.7 

5.6 

7.5 
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TABLE  6a.  Summary  of  results  for  parameter  Lef£  ( Equator  Lai ) • 


Average  ALaff  (*) 


PARAMETER 

-20* 

-13H 

-10* 

-5* 

+  5* 

♦  10* 

♦  15* 

♦  20* 

feF2 

-36.4 

-44.8 

-32.2 

-17.3 

20-0 

43.8 

69.2 

99- 0 

V2 

7.2 

5.1 

3.3 

1.6 

-1.7 

-3.0 

-4.1 

-5.3 

-3.9 

-3.2 

-1.8 

-1.0 

0.8 

1.7 

2.4 

3.0 

Yt 

N.(h.) 

-12.4 

-9.4 

-6.4 

-3.2 

3.3 

6.6 

10.0 

13.5 

51.9 

35.3 

21.7 

10.0 

-8.7 

-16.2 

-22  *  8 

-28.6 

RMS 

(*) 

PARAMETER 

-20* 

•W 

-15* 

-10* 

-3* 

♦  5* 

♦  10* 

♦  15* 

♦  20* 

V2 

2.2 

3.0 

2.2 

1.2 

1.4 

2.8 

4  .6 

6-3 

h«F2 

4.0 

2.6 

1.7 

0.9 

0.7 

1.3 

1.6 

1.9 

h~ 

3.7 

2.7 

2.0 

1.0 

0.9 

2.0 

2.9 

3.8 

Yt 

2.0 

1.3 

1.0 

0-5 

0.5 

1.0 

1.4 

1.9 

N„Ch.) 

3.5 

2.4 

1.5 

o.  / 

0.6 

11 

1.6 

2.0 

Minimum  ALaj£  (*) 


PARAMETER 

-20* 

-13* 

-10* 

-3* 

♦  5* 

♦  10* 

♦  15* 

♦  20* 

V2 

-30.  J 

-36.5 

.26.2 

-14.2 

16.3 

35-  3 

57.1 

82.4 

V2 

4.6 

3.1 

2.1 

0.9 

-1.1 

-1.9 

-2.8 

-3.8 

h” 

-0.3 

-0.2 

-5.7 

00 

0.0 

0  0 

0.0 

0-0 

-0.1 

-6.3 

-4.2 

-2.2 

2.2 

4.5 

6-9 

9.5 

Na(h„) 

42  8 

29.2 

17-7 

8.2 

-7.1 

-13.3 

-18.3 

-23-3 

Maximum  Al#f£  (*) 


PARAMETER 

-20* 

-15* 

-10* 

-5* 

♦  5* 

l  ♦ 

t  >— 

i  o 

1  if 

+  15* 

♦  20* 

V2 

-57.6 

-46.5 

-33.5 

-18.0 

20.8 

44.7 

72.3 

103.5 

hmF2 

17.7 

11.6 

7.5 

3.7 

-3.4 

-5.9 

-8.5 

-10.5 

h? 

-10.6 

-8.4 

0.6 

-3.0 

2.7 

3.7 

8-4 

10-9 

Yt 

-13.1 

-11.5 

-7.7 

-3.9 

3.9 

7.9 

11.9 

16.1 

N,<*V 

54.1 

37.0 

22.7 

10.5 

-9.2 

-16-8 

-23.7 

-29.7 
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TABLE  6b.  Summary  of  raaulta  for  paramatar  <Ltl»j>  (Equatorial) 


Avaraga  A<Lff£>  (X) 


PARAMETER 


f0F2 

6.5 

3.8 

3.3 

1.6 

-1.3 

-2.3 

-3.2 

-A  •  0 

hmF2 

7.2 

3.2 

3. A 

1.7 

-1.6 

-3.0 

-A  •  0 

-3.2 

n! 

-3.9 

-3.1 

-1.8 

-0.9 

0.9 

1.7 

2. A 

3.1 

Yt 

-12.3 

-9.  A 

-6*3 

-3.2 

3.3 

6.6 

10.1 

13.5 

N#(hf) 

-2.7 

-2.0 

-1.4 

-0.7 

0.7 

l.A 

2.2 

2.9 

RMS  A<Lgff> 

(X) 

PARAMETER 

-20X 

-13X 

-10X 

-3X 

♦  5X 

tlOX 

♦  15X 

♦  20X 

fcF2 

3  3 

3.8 

3.3 

15 

1.1 

2.0 

2.6 

3.1 

h°F2 

h? 

AO 

2.6 

1.7 

0.8 

0.8 

l.A 

1.5 

1.9 

3.7 

2.8 

2.0 

1.0 

1.0 

2.0 

2.9 

3.9 

Vt 

2.0 

1.5 

1.0 

0.5 

0.5 

1.0 

1.3 

1.9 

N,<h,) 

2.2 

1.7 

1.2 

0.6 

0.7 

1.3 

2.1 

2.8 

Minimum  A<Lgf£>  (X) 


PARAMETER 

-20X 

-13X 

-10X 

i 

i  ■ 

1  m 

1  * 

♦  3x 

♦  10X 

♦  13x 

♦  20* 

fcF2 

3.7 

2.5 

1.5 

0.7 

-0.6 

-1.1 

-1.5 

-1.9 

A. 6 

3.3 

2.1 

1.0 

-0.9 

-1.9 

-2.7 

-3.6 

h* 

-0.3 

-0.1 

-5.7 

0.0 

0.0 

0.0 

0.0 

0.0 

Yt 

-8.1 

-6.2 

-A. 2 

-2.1 

2.2 

A. 5 

6.9 

9.3 

Ng(hg) 

-1.2 

-0.9 

-0.6 

-0.3 

0.3 

0.6 

1.0 

1.3 

PARAMETER 


W 


Maximum  A<Lr|f>  (X) 


-20X 

■  13X 

-10X 

-5X 

♦  5X 

♦  10X 

♦  15X 

♦  20* 

V 

.  j 

22.3 

21.8 

12.6 

53 

-A.  2 

-7.6 

-10-1 

-12.1 

Si 

17.8 

11.6 

7.7 

3.7 

-3.  A 

-5.9 

-8.3 

-10.5 

H 

-10.6 

-8.3 

0.6 

-2.9 

2.9 

5.7 

8. A 

11.0 

-13.1 

-11. A 

-7.7 

-3.9 

3.9 

8-0 

12.0 

16-1 

|l 

-8-6 

-6.6 

-A. 5 

-2-3 

2.5 

3-1 

7.7 

10.5 

t 

■aft*  rfWadU  afta'.lrfU  ML  mJl-  ±jCmj£>  ±  m  *J(  •  X  WA«,  4~ «.  it  */a  A  *  M  *  .<  #  4.  *4' »  4  •  Jt  -  4  m-A  *,A  *,4  +.A  a»A  *Jk  *«A  •  ",  *.1  *,.1  *J5.  *,1  w\  a^a.'1' 


TABLE  7a.  Summary  of  results  for  parameter  Leff  (All  cases)- 


Average  (*) 


PARAMETER 

-20* 

-15* 

1  I 
|  1— * 

1  o 

1  X 

-5* 

♦  5K 

♦  10* 

♦  13H 

♦  2  0X 

f_P2 

-54. 3 

-43.4 

-31.1 

-16.7 

19.1 

40.9 

65.9 

94.1 

o 

h_F2 

7.0 

5.2 

3.4 

1.6 

-1.7 

-3.2 

-4.5 

-3.8 

h” 

-3.1 

-2.3 

-1.4 

-0.7 

0.6 

1.1 

1.5 

1.9 

nr 

Yt 

-9.3 

-7.1 

-4.8 

-2.4 

2.5 

5.0 

7.7 

10.4 

Ng(ha) 

49.5 

33.9 

20.8 

9.6 

-8.4 

-13.6 

-22.0 

-27.6 

RMS  AL,{{ 

(*> 

PARAMETER 

-20* 

-13* 

-10* 

-3* 

♦  5* 

♦  10* 

♦  13* 

mm  mm  MB 

♦  20* 

f~P2 

2.7 

2.7 

2.0 

1.1 

.  1.3 

2.9 

4.9 

7.1 

0 

h_F2 

3.1 

2.1 

1.3 

0.7 

0.6 

1.1 

1.3 

1.7 

‘  m 

|\ija 

2.3 

2.1 

1.3 

0.8 

0.7 

1.3 

2.2 

2.9 

y| 

3.4 

2.6 

1.7 

0.9 

0.9 

1.8 

2.6 

3.3 

N,(hg) 

3.6 

2.4 

1.3 

0.7 

0.6 

1.0 

1.4 

1.8 

Minimum  (*> 


PARAMETER 

-20* 

-13* 

-10* 

-3* 

♦  5* 

♦  10* 

♦  13* 

♦  20* 

f_F2 

-30.0 

-36.3 

-26.2 

-14.2 

16-3 

33.3 

37.1 

82.4 

o 

h1_F2 

2.4 

2.2 

1.6 

0.9 

-1.0 

-1.9 

-2.8 

-3.8 

Btt 

hT 

-0.3 

-0.2 

-3.7 

o.o 

o.o 

0.0 

0.0 

0.0 

*  T 

Yt 

-4.7 

-3.6 

-2.4 

-1.2 

1.3 

2.3 

4.0 

3  ■  3 

Hg(hg) 

42.8 

29.2 

17.7 

0.2 

-7.1 

-13.3 

-18-5 

—  23*3 

\y 

I'.V 


& 

v> 

fr*‘ 

* 

!?: 

W( 


Mr 

V1 

!»’ 

i1 

V, 

pi:’ 

fe 


I 


Maximum  Al>vff  (X) 


PARAMETER 

-20* 

-13* 

-10* 

-3* 

♦  5* 

♦  10* 

♦  15* 

♦  20* 

f_F2 

-37.6 

-46-3 

-33.3 

-18.0 

20.8 

44.7 

72.3 

103.3 

o 

h_F2 

17.7 

11.6 

7.3 

3.7 

-3.4 

-5.9 

-8.5 

-10.3 

m 

hT 

-10.6 

-8.4 

0.6 

-3.0 

2.7 

5.7 

8.4 

10.9 

yJ 

-15.1 

-11.3 

-7.7 

-39 

3.9 

7.9 

11.9 

16.1 

Nf(h-) 

34.1 

37.0 

22.7 

10.5 

-9.2 

-16.8 

-23-7 

-29-7 

Ki 

5:5 

Sv 


K/ 

£ 

I 


TABLE  7b.  Suaaary  of  reaulta  for  parameter  <L#jj>  (All  casaa). 


Average  A<t.g£f>  (4) 


PARAMETER 

— 20*» 

-134 

-104 

*  i 
«■»  i 
i  i 

♦  34 

♦  104 

♦  134 

♦  204 

V2 

11.6 

8.3 

3.1 

2.3 

-2.0 

-3.7 

-3.1 

-6.4 

V2 

7.1 

3.2 

3.4 

1.7 

-1.6 

-3.1 

-4.4 

-3.7 

hT 

-3.1 

-2  •  3 

-1  >  3 

-0.7 

0.6 

1.1 

1.3 

1.9 

yt 

-9.3 

-7.1 

-4.8 

-2.4 

2.5 

3.1 

7.7 

10-4 

Nfl(hg) 

-4.2 

-3 . 2 

-2 . 1 

-1.1 

1.1 

2.2 

3.2 

4.4 

PARAMETER 

-204 

-134 

RMS 

-104 

*<Leff> 

-34 

(4) 

♦  34 

♦  104 

♦  134 

♦  204 

V2 

6.6 

3.1 

3.0 

1.4 

1.1 

2.0 

2.8 

3.4 

H.P2 

3.1 

2.1 

1.3 

0.6 

0.6 

1.1 

1.3 

1.7 

h- 

2.3 

2.0 

1.3 

0.8 

0.8 

1.3 

2.2 

3.0 

3.4 

2.6 

1.7 

0.9 

0.9 

1.8 

2.6 

3.3 

2.3 

1.8 

1.2 

0.6 

0.6 

1.3 

1.9 

2.6 

Mlnlnua 

A<La£f>  (4) 

PARAMETER 

-204 

-134 

-104 

1 

1  1 

1  Ui 

>  jr 

♦  34 

♦  104 

♦  134 

♦  204 

V2 

3.7 

2.3 

1.3 

0.7 

-0.6 

-1.1 

-1.3 

-1.9 

hBF2 

2.4 

2.2 

1.7 

0.9 

-0.9 

-1.9 

-2.7 

-3.6 

hT 

-0.3 

-0.1 

-3.7 

00 

0.0 

0.0 

0.0 

0.0 

-4.7 

-3.6 

-2.3 

-1.3 

1.3 

2.7 

4.0 

33 

H.CK,) 

-1.2 

-0.9 

-0-6 

-0.3 

0.3 

0.6 

1.0 

1.3 

Maximus 

A<L#ff>  (*) 

PARAMETER 

-204 

-134 

-104 

-54 

♦  54 

♦  104 

♦  154 

♦  204 

'o'2 

22.3 

21.8 

12.6 

33 

-4.2 

-7.6 

-10.1 

-12.1 

hm'2 

17.8 

11.6 

7.7 

3.7 

-3.4 

-3.9 

-8.3 

-10. 3 

h. 

-10.6 

-8.3 

0.6 

-2.9 

2.9 

3.7 

8.4 

11.0 

-13.1 

-11.4 

-7.7 

-3.9 

3.9 

6.0 

12.0 

16.1 

N.(h.) 

-8.6 

-6.6 

-4.3 

-2.3 

2.3 

3.1 

7.7 

10.5 
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to  +  7%.  For  the  Equatorial  cases,  however,  the  ALefj  values  were 
twice  that  of  the  Auroral  cases,  with  a  range  of  -12%  to  +14%. 

(5)  Ne(hs),  The  effects  of  uncertainties  in  this  parameter  on 
I  .  Leff  were  second  only  to  the  effects  of  uncertainties  in  f0F2,  ranging 
|  from  roughly  +50%  to  -30%.  Fortunately,  this  parameter  will  be 
!  .  measured  directly  at  the  satellite  with  a  probable  error  of  +5%  or  so, 

I  which  translates  to  an  uncertainty  in  L#fl»  of  roughly  +10%. 

The  results  for  were  similar,  with  the  notable  exception 
of  the  level  of  effects  due  to  uncertainties  in  the  two  density 
parameters,  f0F2  and  Ne(hB).  With  the  removal  of  the  density  ratio 
from  the  effects  of  +20%  errors  in  these  two  parameters  drop  to 
ranges  of  +12%  to  -6%  and  -4%  to  +4%  for  f0F2  and  Ne(hs), 
respectively.  The  results  for  the  other  parameters  were  identical  to 
those  for  since  both  f0F2  and  NB(h8)  were  held  constant  (or 
nearly  so)  for  these  cases. 

The  main  difficulty  in  using  these  results  is  assessing  what 
level  of  errors  to  expect  in  each  of  the  parameters  and  how  they  may 
be  correlated  with  one  another.  For  example,  errors  in  either  f0F2  or 
NB(h8)  result  in  large  errors  in  both  and  C^L.  However,  if  the 
errors  in  these  parameters  are  strongly  correlated  such  that  the  ratio 
between  Ne(he)  and  NmF2  is  nearly  constant,  the  errors  in  C^L  will  be, 
at  most,  more  on  the  level  of  those  found  in  <Le;tf>.  This  was 
actually  done  for  case  E-07,  i.e.  f0F2  was  stepped  from  -20%  to  +20% 
and  Ne(hs )  was  changed  to  Keep  the  density  ratio  constant,  and  the 
error  in  was  zero.  We  can,  however,  at  least  make  the  following 
observations : 

( 1 )  The  most  crucial  profile  parameters  for  making  accurate 
estimates  of  C^L,  assuming  that  the  variation  of  <ANe2>  with  height 
is  described  by  <ANe2>1/2/NQ  ■  constant,  are  the  electron  density  at 
the  profile  peak  and  at  the  satellite.  Errors  of  +20%  in  these  can 
lead  to  errors  in  C^L  of  a  factor  of  2. 

(2)  Errors  in  hT  will  have  no  appreciable  effect  on  C^L  unless 
the  transition  height  occurs  near  or  below  the  satellite  altitude. 
This  should  only  be  a  problem  during  night  in  the  equatorial  region 
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during  periods  of  low  sunspot  number,  and  the  effect  on  C^L  should  not 
exceed  roughly  10%. 

(3)  The  effect  of  errors  in  the  height  and  shape  of  the  peak 
(hmF2  and  Y*.)  on  Cj<L  will  typically  not  exceed  5-8%. 

(4)  Noting  that  (a)  the  range  of  < I*e c f >  is  also  much  less  than 
that  of  Lgffr  and  (b)  the  lesser  effect  of  errors  in  the  two  density 
parameters  on  <LCff>,  it  may  make  sense  to  develop  a  global  model  for 
<Leff>  and  couple  it  with  observations  of  Ne(he)  and  values  for  fQF2 
from  a  good  model  or  analysis  to  calculate  C^L  rather  than  attempting 
to  model  or  calculate  Ltff  directly. 

In  closing,  it  should  be  remembered  that  these  results  pertain 
only  to  the  case  where  (1)  the  satellite  is  within  the  Irregularity 
layer  and  has  taken  a  sample  representative  of  the  entire  layer,  (2) 
the  height  distribution  of  the  irregularities  can  be  modeled  by 
<dNg2>1/2/Nt  ■  constant  in  a  layer  between  hmF2  and  hB,  and  (3)  the 
height  variation  of  the  electron  density  profile  in  the  topside  can  be 
approximated  by  a  two-component  plasma  (0+  and  H+)  in  diffusive 
equilibrium. 


S.  CONCLUSION 


This  report  presented  the  results  of  studies  conducted  during  the 
first  year  of  this  project  aimed  at  developing  methods  for  calculating 
estimates  of  the  ionospheric  irregularity  parameter  C^L  from  in  situ 
observations  of  the  ionosphere  from  the  DMSP  SSIES  sensor  package.  It 
was  found  that  the  two  methods  for  calculating  C^  from  the  ion  density 
data  had  the  same  levels  of  uncertainty  in  the  final  value  for  as 
long  as  the  data  vers  properly  processed.  The  processing  methods 
found  to  provide  the  best  estimates  for  C^  were  as  follows: 

(1)  T]/q  method.  Detrend  the  data  set  using  (at  least)  a 
quadratic  detrender;  window  the  data  with  a  30%  split-bell  cosine 
taper;  calculate  and  estimate  the  PDS  using  an  PPT  of  the  windowed 
data;  smooth  the  PDS  using  a  5-point,  centered,  binomial-weight 
smoother;  calculate  estimates  of  T^  and  q  from  a  log-linear  least- 
squares  fit  to  the  smoothed  PDS  over  the  frequency  range  0. 5-7.0  Hz. 

(2)  <ANt2>/q  method.  Detrend  the  data  using  a  Fourier-type 
detrender;  calculate  <ANa2>  from  tho  detrended  data  set;  calculate  an 
estimate  for  q  using  the'  detrended  data  and  the  procedures  described 
above  for  the  Tt/q  method. 

The  critical  element  for  the  T^/q  method  is  that  the  data  need  to  be 
windowed  prior  to  the  FFT  and  the  resulting  PDS  smoothed  prior  to  the 
log-linear  fit  to  obtain  the  best  estimates  for  T1  and  q.  For  the 
<ANa2>/q  method,  the  critical  elements  are  to  (1)  use  a  frequency- 
domain  filter  for  detrending  so  that  the  cutoff  frequency  is  well 
defined,  and  (2)  use  an  observed  value  for  q  rather  than  a  mean,  or 
modeled,  value.  Using  these  processing  methods,  the  resulting 
uncertainties  in  due  to  uncertainties  in  parameters  obtained  from 
the  density  data  set  are  on  the  order  of  5  to  10%. 

The  effects  of  uncertainties  in  the  calculated  value  of  the 
effective  probe  velocity,  Vp,  on  Cj<  were  also  investigated,  ana  were 
found  to  produce  uncertainties  in  in  a  rang^  from  a  few  percent  to 
over  100%.  The  main  sources  of  uncertainty  in  Vp  are  the  shape  of  the 
irregularities  (rod-like  or  sheet-like)  and  the  in  zitu  drift  velocity 
of  the  irregularities. 


The  total  level  of  uncertainties  in  from  all  sources  was  found 
to  be  as  follows: 

(1)  Equatorial.  Uncertainties  should  be  on  the  order  of  5  to  10% 
when  the  in  aitu  drift  velocities  are  known,  and  10-20%  when  they  are 
not. 

(2)  Auroral/polar .  Uncertainties  should  be  on  the  order  of  5  to 
10%  when  the  in  situ  drift  velocities  and  the  irregularity  shape  are 
well  known,  IS  to  30%  when  the  drift  velocities  are  not  well  known, 
and  25  to  >100%  when  neither  are  well  known. 

In  both  regimes,  the  uncertainty  in  will  increase  roughly  linearly 
with  increasing  uncertainty  in  vp.  The  largest  source  of  uncertainty 
is  the  shape  of  the  irregularities  at  auroral/polar  latitudes.  Since 
this  cannot  be  directly  measured  from  the  SSZES  data  set,  a  good  model 
must  be  provided  for  the  two  axial  ratio  parameters,  a  and  b. 

A  method  was  developed  for  calculating  an  estimate  for  C^L  from 
Cfc  at  the  satellite  altitude.  Two  "layer  thickness"  parameters  were 
defined,  an  effectlvo  layer  thickness,  Ltff,  defined  such  that 

ckL  ■  ckxLeff 

and  a  normalized  effective  layer  thickness,  <Ltfj>>,  defined  by 

<Lef f  >  "  £Ne*hs  )/NmF2]*Laf f • 

Both  axe  functions  of  the  altitude  range  of  the  irregularity  layer  and 
the  altitude  variation  of  <dNa2>  within  the  layer.  A  parametric 
study  or  uncertainties  in  and  <Leff>  was  conducted  in  which  the 
irregularity  layer  was  assumed  to  extend  from  the  F2-layer  peak  to  the 
altitude  of  the  satellite,  <ANe2>1 /2/Ne  was  assumed  to  be  constant 
with  altitude  throughout  the  layer,  and  the  background  ionospheric 
electron  density  was  modeled  by  an  adjustable  diffusive-equilibrium 
distribution. 


The  major  results  of  this  study  were: 

( 1  )  The  most  crucial  profile  parameters  for  making  accurate 
estimates  of  CkL  are  the  electron  density  at  the  F2-layer  peak  and  at 
the  satellite  altitude.  Errors  of  only  20%  in  f0F2  can  result  in  100% 
errors  in  CkL. 

(2)  The  effects  of  errors  in  the  densities  are  much  less  on 
<Leff>  du«  to  the  decoupling  of  the  densities  from  this  parameter • 

(3)  The  effects  of  +20%  errors  in  all  other  profile  parameters 
resulted  in  less  than  10%  errors  in  CkL. 

Based  on  this  study,  future  development  of  techniques  for  converting 
Ck  to  CkL  will  focus  on  developing  a  global  model  for  <L.ejf>  which  can 
be  coupled  to  an  external  source  for  the  densities  to  calculate 
This  investigation  of  the  uncertainties  in  the  conversion  of  Ck  to  ckL 
will  be  continued  in  the  next  year  with  a  study  of  the  effects  of 
uncertainties  in  the  height  distribution  of  the  irregularities. 
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Appendix  A.  Data  Sets  Used  in  Studies 


A.1  Simulated  Data  Sets 

In  order  to  measure  the  effectiveness  and  accuracy  of  the  various 
techniques  used  in  calculating  from  a  plasma  density  data  sample, 
two-  sets  of  simulated  density  data  with  known  spectral  characteristics 
were  generated.  Such  a  data  set  can  be  constructed  simply  from 
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aj_  sin 
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[A- 1  ] 


where  ANn  is  the  nth  density  data  point  in  the  sample,  N  is  the  total 
number  of  points  per  data  sample,  Ai  is  a  frequency  shift  parameter, 
is  a  random  phase  shift  parameter  (range:  -ir  to  +*),  and  a^  is  the 
amplitude.  The  frequency  shift  parameter,  Ai,  subtracted  from  i  in 
calculating  the  argument  of  the  sine  function,  is  used  to  shift  the 
power  in  the  data  set  to  frequencies  away  from  frequencies  of  the 
discrete  FFT  used  to  calculate  estimates  of  the  spectral  density 
functions  from  the  data  sets.  This  is  done  to  maximize  the  effects  of 
spectral  "leakage"  between  spectral  frequency  bins.  For  example  if 
Ai  ■  0,  all  power  is  located  at  the  FFT  frequencies  and  no  leakage 
will  occur.  If,  on  the  other  hand,  Ai  •  1/2,  all  power  is  located  at 
frequencies  halfway  between  the  FFT  frequencies,  which  will  provide  a 
measure  of  the  maximum  spectral  leakage. 

The  amplitude,  a^,  is  given  by 

ai  “  LA-2 J 

where  Af  is  the  frequency  step  ( 2  4  /  N  for  SSIES  data  sets),  and  is 
the  desired  power  density  spectrum  ( PDS )  for  the  sample.  For  our 
studies,  we  assume  that  the  PDS  can  be  described  by  a  power  law  of  the 
form 


where  T1  is  the  PDS  value  at  a  frequency  of  1  Hz  and  q  is  the  slope  of 
the  PDS.  The  frequencies  to  use  in  Equation  [A-3]  are  given  by 
m  (i-Ai)Af,  so  combining  Equations  [A-2]  and  [A-3],  the  amplitude 
then  becomes 


ai  ■  2N/T,  (i-Ai)-<3. 


[A-4] 


Ae  we  may  want  to  set  up  different  samples  for  different  levels 
of  irregularity  strength  (T1  ),  we  define  a  normalized  density  sample 
by 


V  ' 


<ANn>  « 


[A-5J 


Using  Equations  [A-l],  [A-4],  and  [A-5],  the  normalized  density  is 

calculated  from 


<ANn>  *  2Af“^-1  )/2  jr  ( i-Ai  )-<?/2sin  (i-Ai)-jj-  +1^  . 

i-1 
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We  can  also  define  a  normalized  RMSAn,,,  <RMS>,  from 


<  RMS  > 


£ANn2W2 


«  T,1/2  RMSANn 


[A-7  ] 


Therefore,  given  q,  we  can  specify  a  normalized  density  sample  and 
RMS;  and  given  T ^ ,  we  can  calculate  the  desired  density  sample  and  the 
RMSANn  for  the  sample  from 


ANn  ■  T, 1 /2  < ANn  > 


[A-8] 


[A-9] 


RMSAN  =  1 <RMS>  . 

In  the  studies  presented  here,  two  sample  data  bases  were 
constructed  containing  six  2048-point  realizations  for  each  q  value  in 
the  set  q:q  *  1.0,  1.2,  1.4,  ...,  2.8,  3.0  ,  one  data  base  for  the 
zero-leakage  case,  (Ai  ■  1/2).  Figure  A-1  shows  one  of  the 
realizations  from  each  data  base  for  q  -  1.6,  and  Figure  A-2  shows  the 
PDS  constructed  from  each  realization  (note  that  the  PDS  for  the  zero- 
leakage  case  is  shifted  up  a  decade  on  the  plot).  No  windowing  or 
detrending  was  done  in  constructing  these  PDS  estimates,  and  the 
effects  of  the  spectral  leakage  from  the  maximum-leakage  case  are 
quite  evident.  The  variation  of  RMS  as  a  function  of  q  for  the  two 
simulation  data  bases  is  shown  in  Figure  A-3. 

In  order  to  study  detrender  effects,  it  was  decided  to  add  more 
terms  to  Equation  [A -6]  with  wavelengths  longer  than  the  data  sample 
size.  The  frequencies  for  these  additional  terms  were  selected  as 
0 . 25Af ,  0 . 5Af ,  and  075A1.  All  three  are  added  to  the  zero-leakage 
cases,  but  only  the  0.25Af  and  0.75Af  terms  are  added  to  the  maximum- 
leakage  cases,  as  the  0.5Af  term  has  already  been  included  in 
Equation  [A-6] .  These  terms  were  added  on  as  density  trends  during 
the  parametric  studies  rather  than  directly  to  the  simulation  data 
bases.  This  was  done  so  that  the  zero-leakage  data  samples  could  be 
used  as  "sanity  checks"  for  the  processing  software,  since  they  should 
return  a  straight-line  PDS  when  the  data  are  neither  detrended  nor 
windowed,  as  in  Figure  A-2. 

A. 2  WIDEBAND  Phase  Scintillation  Data  Set 

No  in  situ  plasma  density  data  sets  were  available  for  this 
study.  Fortunately,  a  large  body  of  phase  scintillation  data  is 
available  at  NWP.A  from  both  the  WIDEBAND  and  HiLat  beacon  experiments. 
While  these  data  are  not  direct  measures  of  the  one-dimensional  in 
situ  density  spectrum,  they  do  provide  a  measure  of.  the  two- 
dimensiona.l  spectrum  and  have  similar  characteristics,  i.e.,  they  can 
be  approximated  by  a  red-noise  power-law  spectrum. 
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Figure  A-2.  Power  density  spectra  of  the  data  In  Figure  A-l. 
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Figure  A-3.  Variation  of  RMS  with  q  for  the  zero-  and  maximum- 
leakage  data  sets. 
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A  single  WIDEBAND  satellite  pass  at  the  Poker  Flat,  Alaska, 
receiving  site  was  selected  for  this  study  (pass  PF-52-47).  Figure  A- 
4  shows  plots  of  the  phase  data  from  this  pass.  The  upper  plot  in 
this  figure  is  the  raw  phase  record  for  the  pass.  The  lower  plot  is  a 
detrended  phase  obtained  by  removing  a  low-pass  trend  calculated  using 
a  6-pole  Butterworth  filter  similar  to  one  used  in  routine  WIDEBAND 
processing.  In  this  case,  the  6dB  cutoff  frequency  was  set  at 
0.03333Hz  (a  30  second  detrend  time).  For  this  study,  the  phase  data 
are  reduced  from  the  full  sampling  rate  of  500Hz  down  to  25Hz.  A 
total  of  64  data  samples  are  extracted  from  the  pass  by  selecting  512- 
point  data  samples  at  256-point  intervals  throughout  the  pass. 


(RAD)  PHASE  (RAO) 


POKER  FLAT  PA88  PF-52-47 


180 


100 


1225 


1227 


1228 


1231 


1233 


GMT 


-J- 


1235 


1237 


Figure  A-4.  Phase  scintillation  from  WIDEBAND/Poker  Flat  pass  PF-52-47. 

Upper  plot  Is  the  raw  phase  record;  lower  plot  Is  the  detrended 
phase  obtained  by  removing  a  low-pass  filter  trend. 
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Appendix  B.  Topside  Electron  Density  Model 


The  electron  density  profile  model  used  in  this  study  was 
developed  specifically  fox  use  with  the  DMSP  8S1ES  sensor  package. M 
The  underlying  structure  of  the  topside  section  of  the  model  is  that 
of  a  two-component  ionosphere  (04  and  H4)  in  diffusive  equilibrium, 
but  it  has  been  parameterized  to  allow  the  model  to  be  fit  to  non¬ 
equilibrium  conditions.  Assuming  charge  neutrality,  the  electron 
density  profile  will  be  identical  to  the  height  variation  of  the 
ionospheric  plasma  density,  which  for  this  model  is  given  by 


N#(h)  ■  Np(h)  « 

a  ■  ckq  4  a1  (h-400 ) 


rN 

\  Tp  /  ' 


o(04)a"l6rtI  4  H0(H4)e-Ij 


[B-1] 

[B-2] 


where  Tp  is  the  plasma  temperature  (T^4T# )  ;  Tp0  is  the  plasma 


temperature  at  a  reference  height  (h0);  No(04)  and  N0(H4)  are  the 
number  densities  of  04  and  H4  at  the  roference  height;  ag*  and  0 
are  profile  adjustment  parameters;  and  n  and  I  are  integral  functions 
given  by 


/u(  h ) 


dh 


[B-3] 


1(h) 


h  m  „ . 
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dh 
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where  m.  is  the  mass  of  H4,  and  R  is  defined  as  the  ratio  of  N(04)  to 


N(H4) .  Ion  and  electron  temperatures  are  obtained  from  a  lecent  model 
based  on  an  analysis  of  data  from  the  AE-C  satellite .W  The  density 
ratio,  R,  is  based  on  the  04  to  H4  transition  height,  hr,  defined  as 
the  height  at  which  N(04)  ■  N(H+'.  This  parameter  can  be  either  (1) 
calculated  from  observations  of  N(04),  N(H4),  and  Tp  from  the  SSIES 


B-1 


RPA  sensor,  or  (2)  obtained  from  an  empirical  model  of  hT  derived  from 
published  analyses  of  topside  profiles  from  the  Alouette  satellites 
and  RPA  data  from  0G0-6.1®! 


The  topside  model  is  fit  to  the  F2  peak  by  means  of  a  parabolic 
layer  tAken  from  the  Bent  profile  modelt10)  of  the  form 


Nt(h)  ■  NmF2 


[B-5] 


where  NmF2  and  hmF2  are  the  density  and  height  of  the  F2  layer  peak 
and  Yt  is  the  parabolic  semi-thickness.  The  Y^.  parameter  can  be 
either  estimated  in  the  procedure  used  to  fit  the  profile  to 
observations,  or  obtained  from  the  expressions  used  in  the  original 
Bent  model.  Equations  [B-l]  and  [B-5]  are  fit  together  at  the  height 
where  the  plasma  scale  height  calculated  from  the  two  representations 
of  Np(h)  are  equal.  This  must  be  calculated  iteratively,  but  rarely 
requires  more  that  four  or  five  iterations.  This  height  is  also  used 
as  the  reference  height,  h0,  for  the  parameters  in  Equation  [B-l]. 

For  the  sake  of  providing  a  complete  plasma  density  profile,  the 
bottomside  section  of  the  Air  Weather  Service  (AWS)  RBTEC  model  is 
used  to  describe  the  height  variation  below  the  F2  peak.l,1,2l  This 
model  uses  three  Chapman-function  layers  to  describe  the  three  main 
ionospheric  layers  (E,  FI,  and  F2 ) .  The  choice  of  models  for  the 
bottomside  of  the  profile  has  little  impact  on  the  present  study,  as 
the  irregularity  layer  is  assumed  start  either  at  or  just  slightly 
below  the  F2  peak. 

A  number  of  techniques  were  devised  for  fitting  this  profile 
model  to  a  wide  range  of  input  data.  For  the  application  at  hand,  we 
will  assume  that  an  observation  of  Np  at  the  satellite  altitude  is 
available,  and  that  an  observation  of  f0F2  and  hmF2  may  or  may  not  be 
available.  If  all  three  are  available,  the  profile  can  be  adjusted  to 
fit  all  three  by  adjusting  any  two  of  the  four  parameters,  Yt,  ii ,  n0, 
or  ■«].  If  either  f0F2  or  hmF2  is  missing,  a  profile  can  be  fit  to  the 


density  at  the  satellite  and  the  available  F2-layer  parameter 
providing  an  estimate  of  the  missing  F2-layer  parameter  assuming  fixed 
values  for  Y^  and  the  three  adjustment  parameters. 
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